False gene and chromosome losses in genome assemblies caused by GC content variation and repeats.

BACKGROUND: Many short-read genome assemblies have been found to be incomplete and contain mis-assemblies. The Vertebrate Genomes Project has been producing new reference genome assemblies with an emphasis on being as complete and error-free as possible, which requires utilizing long reads, long-range scaffolding data, new assembly algorithms, and manual curation. A more thorough evaluation of the recent references relative to prior assemblies can provide a detailed overview of the types and magnitude of improvements. RESULTS: Here we evaluate new vertebrate genome references relative to the previous assemblies for the same species and, in two cases, the same individuals, including a mammal (platypus), two birds (zebra finch, Anna's hummingbird), and a fish (climbing perch). We find that up to 11% of genomic sequence is entirely missing in the previous assemblies. In the Vertebrate Genomes Project zebra finch assembly, we identify eight new GC- and repeat-rich micro-chromosomes with high gene density. The impact of missing sequences is biased towards GC-rich 5'-proximal promoters and 5' exon regions of protein-coding genes and long non-coding RNAs. Between 26 and 60% of genes include structural or sequence errors that could lead to misunderstanding of their function when using the previous genome assemblies. CONCLUSIONS: Our findings reveal novel regulatory landscapes and protein coding sequences that have been greatly underestimated in previous assemblies and are now present in the Vertebrate Genomes Project reference genomes.

Characterization of the complete chloroplast genome of Brassica oleracea var. italica and phylogenetic relationships in Brassicaceae.

Broccoli (Brassica oleracea var. italica) is an important B. oleracea cultivar, with high economic and agronomic value. However, comparative genome analyses are still needed to clarify variation among cultivars and phylogenetic relationships within the family Brassicaceae. Herein, the complete chloroplast (cp) genome of broccoli was generated by Illumina sequencing platform to provide basic information for genetic studies and to establish phylogenetic relationships within Brassicaceae. The whole genome was 153,364 bp, including two inverted repeat (IR) regions of 26,197 bp each, separated by a small single copy (SSC) region of 17,834 bp and a large single copy (LSC) region of 83,136 bp. The total GC content of the entire chloroplast genome accounts for 36%, while the GC content in each region of SSC,LSC, and IR accounts for 29.1%, 34.15% and 42.35%, respectively. The genome harbored 133 genes, including 88 protein-coding genes, 37 tRNAs, and 8 rRNAs, with 17 duplicates in IRs. The most abundant amino acid was leucine and the least abundant was cysteine. Codon usage analyses revealed a bias for A/T-ending codons. A total of 35 repeat sequences and 92 simple sequence repeats were detected, and the SC-IR boundary regions were variable between the seven cp genomes. A phylogenetic analysis suggested that broccoli is closely related to Brassica oleracea var. italica MH388764.1, Brassica oleracea var. italica MH388765.1, and Brassica oleracea NC_0441167.1. Our results are expected to be useful for further species identification, population genetics analyses, and biological research on broccoli.

Genome-Wide Prediction of Transcription Start Sites in Conifers.

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

Human SARS-CoV-2 has evolved to increase U content and reduce genome size.

Infections caused by SARS-CoV-2 have brought great harm to human health. After transmission for over two years, SARS-CoV-2 has diverged greatly and formed dozens of different lineages. Understanding the trend of its genome evolution could help foresee difficulties in controlling transmission of the virus. In this study, we conducted an extensive monthly survey and in-depth analysis on variations of nucleotide, amino acid and codon numbers in 311,260 virus samples collected till January 2022. The results demonstrate that the evolution of SARS-CoV-2 is toward increasing U-content and reducing genome-size. C, G and A to U mutations have all contributed to this U-content increase. Mutations of C, G and A at codon position 1, 2 or 3 have no significant difference in most SARS-CoV-2 lineages. Current viruses are more cryptic and more efficient in replication, and are thus less virulent yet more infectious. Delta and Omicron variants have high mutability over other lineages, bringing new threat to human health. This trend of genome evolution may provide a clue for tracing the origin of SARS-CoV-2, because ancestral viruses should have lower U-content and probably bigger genome-size.

Complete Chloroplast Genome Sequence of Sonchus brachyotus Helps to Elucidate Evolutionary Relationships with Related Species of Asteraceae.

Sonchus brachyotus DC. possesses both edible and medicinal properties and is widely distributed throughout China. In this study, the complete cp genome of S. brachyotus was sequenced and assembled. The total length of the complete S. brachyotus cp genome was 151,977 bp, including an LSC region of 84,553 bp, SSC region of 18,138 bp, and IR region of 24,643 bp. Sequence analyses revealed that the cp genome encoded 132 genes, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The GC content was 37.6%. One hundred mononucleotide microsatellites, 4 dinucleotide microsatellites, 67 trinucleotide microsatellites, 4 tetranucleotide microsatellites, and 1 long repeat were identified. The SSR frequency of the LSC region was significantly greater than that of the IR and SSC regions. In total, 175 SSRs and highly variable regions were recognized as potential cp markers. By analyzing the IR/LSC and IR/SSC boundaries, structural differences between S. brachyotus and 6 other species were detected. According to phylogenetic analyses, S. brachyotus was most closely related to S. arvensis and S. oleraceus. Overall, this study provides complete cp genome resources for S. brachyotus that will be beneficial for identifying potential molecular markers and evolutionary patterns of S. brachyotus and its closely related species.

The Second Chromosome Promotes the Adaptation of the Genus Flammeovirga to Complex Environments.

Approximately 10% of bacterial strains contain more than one chromosome; however, in contrast to the primary chromosomes, the mechanisms underlying the formation of the second chromosomes and the significance of their existence remain unclear. Species of the genus Flammeovirga are typical polysaccharide-degrading bacteria, and herein, we report complete genome maps of this genus. These genomes all had multireplicons and second chromosomes. The second chromosome, much larger than plasmids and even megaplasmids, had rRNA and a disparity of 1% relative to the main chromosome in guanine-cytosine (GC) content. The largest chromosomes carried core genes for cellular processes, while the second chromosomes were enriched with genes involved in the transport and metabolism of inorganic ions and carbohydrates, particularly genes encoding glycoside hydrolases and polysaccharide lyases, which constituted the genetic basis for the strains' excellent capabilities to utilize polysaccharides. The second chromosomal evolution had a higher mutation rate than the primary chromosomes. Furthermore, the second chromosomes were also enriched in horizontal transfer genes and duplicated genes. The primary chromosomes were more evolutionarily conserved, while the second chromosomes were more plastic, which might be related to their different roles in the bacterial survival process. This study can be used as an example to explain possible formation mechanisms and functions of the second chromosomes, providing a reference for peer research on the second chromosomes. In particular, the second chromosomes were enriched in polysaccharide-degrading enzymes, which will provide theoretical support for using genomic data to mine tool-type carbohydrase resources. IMPORTANCE For decades, the typical bacterial genome has been thought to contain a single chromosome and a few small plasmids carrying nonessential genes. However, an increasing number of secondary chromosomes have been identified in various bacteria (e.g., plant symbiotic bacteria and human pathogens). This study reported three complete genomes of the polysaccharide-degrading marine bacterial genus Flammeovirga, revealed that they harbor two chromosomes, and further identified that the presence of a multireplicon system is a characteristic of complete Flammeovirga genomes. These sequences will add to our knowledge on secondary chromosomes, especially within Bacteroidetes. This study indicated that the second chromosomes of the genus Flammeovirga initially originated from an ancestral plasmid and subsequently expanded by gene duplication or by obtaining heterologous genes with functions, thus promoting host strains to adapt to complex living environments (e.g., to degrade more diverse polysaccharides from marine environments). These findings will promote the understanding of the evolution and function of bacteria with multireplicon systems.

The First Two Complete Mitochondrial Genomes of Neoephemeridae (Ephemeroptera): Comparative Analysis and Phylogenetic Implication for Furcatergalia.

Mayflies of the family Neoephemeridae are widespread in the Holarctic and Oriental regions, and its phylogenetic position is still unstable in the group Furcatergalia (mayflies with fringed gills). In the present study, we determined the complete mitogenomes of two species, namely Potamanthellus edmundsi and Pulchephemera projecta, of this family. The lengths of two mitogenomes were 15,274 bp and 16,031 bp with an A + T content of 73.38% and 73.07%, respectively. Two neoephemerid mitogenomes had a similar gene size, base composition, and codon usage of protein-coding genes (PCGs), and the sequenced gene arrangements were consistent with the putative ancestral insect mitogenomes as understood today. The most variable gene of Furcatergalia mitogenomes was ND2, while the most conserved gene was COI. Meanwhile, the analysis of selection pressures showed that ND6 and ATP8 exhibited a relaxed purifying selection, and COI was under the strongest purifying selection. Phylogenetic trees reconstructed based on two concatenated nucleotide datasets using both maximum likelihood (ML) and Bayesian inference (BI) estimations yielded robust identical topologies. These results corroborated the monophyly of seven studied families and supported the family Leptophlebiidae as being of the basal lineage of Furcatergalia. Additionally, the sister-group relationship of Caenidae and Neoephemeridae was well supported. Methodologically, our present study provides a general reference for future phylogenetic studies of Ephemeroptera at the mitogenome level.

First Description of the Mitogenome and Phylogeny of Culicinae Species from the Amazon Region.

The Culicidae family is distributed worldwide and comprises about 3587 species subdivided into the subfamilies Anophelinae and Culicinae. This is the first description of complete mitochondrial DNA sequences from Aedes fluviatilis, Aedeomyia squamipennis, Coquillettidia nigricans, Psorophora albipes, and Psorophora ferox. The mitogenomes showed an average length of 15,046 pb and 78.02% AT content, comprising 37 functional subunits (13 protein coding genes, 22 tRNAs, and two rRNAs). The most common start codons were ATT/ATG, and TAA was the stop codon for all PCGs. The tRNAs had the typical leaf clover structure, except tRNASer1. Phylogeny was inferred by analyzing the 13 PCGs concatenated nucleotide sequences of 48 mitogenomes. Maximum likelihood and Bayesian inference analysis placed Ps. albipes and Ps. ferox in the Janthinosoma group, like the accepted classification of Psorophora genus. Ae. fluviatilis was placed in the Aedini tribe, but was revealed to be more related to the Haemagogus genus, a result that may have been hampered by the poor sampling of Aedes sequences. Cq. nigricans clustered with Cq. chrysonotum, both related to Mansonia. Ae. squamipennis was placed as the most external lineage of the Culicinae subfamily. The yielded topology supports the concept of monophyly of all groups and ratifies the current taxonomic classification.

Phylogenetic Distribution of WhiB- and Lsr2-Type Regulators in Actinobacteriophage Genomes.

Viruses that infect different actinobacterial host species are known as actinobacteriophages. They are composed of highly divergent and mosaic genomes due to frequent gene exchange between their bacterial hosts and related viral species. This is also reflected by the adaptive incorporation of host transcription factors (TFs) into phage regulatory networks. Previous studies discovered Lsr2-type and WhiB-type regulators encoded by actinobacteriophage genomes. However, limited information is available about their distribution, evolution, and impact on host species. In this study, we computationally screened the distribution of known bacterial and phage TFs inside 2951 complete actinobacteriophage genomes and identified 13 different TF domains. Among those, WhiB, Lsr2, MerR, and Cro/CI-like proteins were widespread and found in more than 10% of the analyzed actinobacteriophage genomes. Neighboring genomic context analysis of the whiB and lsr2 loci showed group-specific conservation of gene synteny and potential involvement of these genes in diverse regulatory functions. Both genes were significantly enriched in temperate phages, and the Lsr2-encoding genomes featured an overall lower GC content. Phylogenetic analysis of WhiB and Lsr2 proteins showed the grouping of phage sequences within bacterial clades, suggesting gene acquisition by phages from their bacterial host species or by multiple, independent acquisition events. Overall, our study reports the global distribution of actinobacteriophage regulatory proteins and sheds light on their origin and evolution. IMPORTANCE Actinobacteriophages are viruses that infect bacterial species of the diverse phylum of Actinobacteria. Phages engage in a close relationship with their bacterial host. This is also reflected by the adoption of genetic material from their host and its incorporation into phage regulatory circuits. In this study, we systematically searched the genomes of actinobacteriophages for the presence of transcription factor domains. We show that proteins belonging to the regulator families of WhiB and Lsr2 belong to the most abundant regulatory proteins encoded by actinobacteriophages. Further phylogenetic analysis shed light on their origin and evolution. Altogether, this study provides an important basis for further experimental investigation of their role in the coordination of the phage life cycle and their interaction with the host regulatory network in this important bacterial phylum.

High-quality whole-genome sequence analysis of Lactobacillus paragasseri UBLG-36 reveals oxalate-degrading potential of the strain.

Lactobacillus paragasseri was identified as a novel sister taxon of L. gasseri in 2018. Since the reclassification of L. paragasseri, there has been hardly any report describing the probiotic properties of this species. In this study, an L. paragasseri strain UBLG-36 was sequenced and analyzed to determine the molecular basis that may confer the bacteria with probiotic potential. UBLG-36 was previously documented as an L. gasseri strain. Average nucleotide identity and phylogenomic analysis allowed accurate taxonomic identification of UBLG-36 as an L. paragasseri strain. Analysis of the draft genome (~1.94 Mb) showed that UBLG-36 contains 5 contigs with an average G+C content of 34.85%. Genes essential for the biosynthesis of bacteriocins, adhesion to host epithelium, stress resistance, host immunomodulation, defense, and carbohydrate metabolism were identified in the genome. Interestingly, L. paragasseri UBLG-36 also harbored genes that code for enzymes involved in oxalate catabolism, such as formyl coenzyme A transferase (frc) and oxalyl coenzyme A decarboxylase (oxc). In vitro oxalate degradation assay showed that UBLG-36 is highly effective in degrading oxalate (averaging more than 45% degradation), a feature that has not been reported before. As a recently identified bacterium, there are limited genomic reports on L. paragasseri, and our draft genome sequence analysis is the first to describe and emphasize the probiotic potential and oxalate degrading ability of this species. With results supporting the probiotic functionalities and oxalate catabolism of UBLG-36, we propose that this strain is likely to have immense biotechnological applications upon appropriate characterization.

Distribution, expression and methylation analysis of positively selected genes provides insights into the evolution in Brassica rapa.

It is believed that positive selection is one of the major evolutionary forces underlying organism phenotypic diversification. Nevertheless, the characteristics of positively selected genes (PSGs), have not been well investigated. In this study, we performed a genome-wide analysis of orthologous genes between Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea), and identified 468 putative PSGs. Our data show that, (1) PSGs are enriched in plant hormone signal transduction pathway and the transcription factor family; (2) PSGs are significantly lower expressed than randomly selected non-PSGs; (3) PSGs with tissue specificity are significantly higher expressed in the callus and reproductive tissues (flower and silique) than in vegetable tissues (root, stem and leaf); (4) the proportion of PSGs is positively correlated with the number of retained triplication gene copies, but the expression level of PSGs decay with the increasing of triplication gene copies; (5) the CG and CHG methylation levels of PSGs are significantly higher in introns and UTRs than in the promoter and exon regions; (6) the percent of transposable element is in proportion to the methylation level, and DNA methylation (especially in the CG content) has the tendency to reduce the expression of PSGs. This study provides insights into the characteristics, evolution, function, expression and methylation of PSGs in B. rapa.

Analysis of Codon Usage Patterns of Six Sequenced Brachypodium distachyon Lines Reveals a Declining CG Skew of the CDSs from the 5'-ends to the 3'-ends.

Brachypodium distachyon, a new monocotyledonous model plant, has received wide attention in biological research due to its small genome and numerous genetic resources. Codon usage bias is an important feature of genes and genomes, and it can be used in transgenic and evolutionary studies. In this study, the nucleotide compositions and patterns of codon usage bias were calculated using Codon W. Additionally, an ENC plot, Parity rule 2 and correspondence analyses were used to explore the major factors influencing codon usage bias patterns. The numbers of hydrogen bonds and skews were used to analyze the GC trend in the 5'-ends of the coding sequences. The results showed that minor differences in the codon usage bias patterns were revealed by the ENC plot, Parity rule 2 and correspondence analyses. The analyses of the CG-skew and the number of hydrogen bonds showed a declining trend in the number of cytosines at the 5'-ends of the CDSs (from the 5'-ends to the 3'-ends), indicating that GC may play a major role in codon usage bias. In addition, our results laid a foundation for the study of codon usage bias patterns in Brachypodium genus and suggested that the GC plays a major role in determining these patterns.

Insights into the Lignocellulose-Degrading Enzyme System of Humicola grisea var. thermoidea Based on Genome and Transcriptome Analysis.

Humicola grisea var. thermoidea is a thermophilic ascomycete and important enzyme producer that has an efficient enzymatic system with a broad spectrum of thermostable carbohydrate-active (CAZy) enzymes. These enzymes can be employed in lignocellulose biomass deconstruction and other industrial applications. In this work, the genome of H. grisea var. thermoidea was sequenced. The acquired sequence reads were assembled into a total length of 28.75 Mbp. Genome features correlate with what was expected for thermophilic Sordariomycetes. The transcriptomic data showed that sugarcane bagasse significantly upregulated genes related to primary metabolism and polysaccharide deconstruction, especially hydrolases, at both pH 5 and pH 8. However, a number of exclusive and shared genes between the pH values were found, especially at pH 8. H. grisea expresses an average of 211 CAZy enzymes (CAZymes), which are capable of acting in different substrates. The top upregulated genes at both pH values represent CAZyme-encoding genes from different classes, including acetylxylan esterase, endo-1,4-ß-mannosidase, exoglucanase, and endoglucanase genes. For the first time, the arsenal that the thermophilic fungus H. grisea var. thermoidea possesses to degrade the lignocellulosic biomass is shown. Carbon source and pH are of pivotal importance in regulating gene expression in this organism, and alkaline pH is a key regulatory factor for sugarcane bagasse hydrolysis. This work paves the way for the genetic manipulation and robust biotechnological applications of this fungus. IMPORTANCE Most studies regarding the use of fungi as enzyme producers for biomass deconstruction have focused on mesophile species, whereas the potential of thermophiles has been evaluated less. This study revealed, through genome and transcriptome analyses, the genetic repertoire of the biotechnological relevant thermophile fungus Humicola grisea. Comparative genomics helped us to further understand the biology and biotechnological potential of H. grisea. The results demonstrate that this fungus possesses an arsenal of carbohydrate-active (CAZy) enzymes to degrade the lignocellulosic biomass. Indeed, it expresses more than 200 genes encoding CAZy enzymes when cultivated in sugarcane bagasse. Carbon source and pH are key factors for regulating the gene expression in this organism. This work shows, for the first time, the great potential of H. grisea as an enzyme producer and a gene donor for biotechnological applications and provides the base for the genetic manipulation and robust biotechnological applications of this fungus.

Complete mitochondrial genome of Dacus vijaysegarani and phylogenetic relationships with congeners and other tephritid fruit flies (Insecta: Diptera).

BACKGROUND: Tephritid fruit flies of the genus Dacus are members of the tribe Dacini, subfamily Dacinae. There are some 274 species worldwide, distributed in Africa and the Asia-Pacific. To date, only five complete mitochondrial genomes (mitogenomes) of Dacus fruit flies have been published and are available in the GenBank. METHODS AND RESULTS: In view of the lack of study on their mitogenome, we sequenced (by next generation sequencing) and annotated the complete mitogenome of D. vijaysegarani from Malaysia to determine its features and phylogenetic relationship. The whole mitogenome of D. vijaysegarani has identical gene order with the published mitogenomes of the genus Dacus, with 13 protein-coding genes, two rRNA genes, 22 tRNAs, a non-coding A + T rich control region, and intergenic spacer and overlap sequences. Phylogenetic analysis based on 15 mitochondrial genes (13 PCGs and two rRNA genes), reveals Dacus, Zeugodacus and Bactrocera forming a distinct clade. The genus Dacus forms a monophyletic group in the subclade containing also the Zeugodacus group; this Dacus-Zeugodacus subclade is distinct from the Bactrocera subclade. D. (Mellesis) vijaysegarani forms a lineage with D. (Mellesis) trimacula in the subcluster containing also the lineage of D. (Mellesis) conopsoides and D. (Callantra) longicornis. D. (Dacus) bivittatus and D. (Didacus) ciliatus form a distinct subcluster. Based on cox1 sequences, the Malaysia and Vietnam taxa of D. vijaysegarani may not be conspecific. CONCLUSIONS: Overall, the mitochondrial genome of D. vijaysegarani provided essential molecular data that could be useful for further studies for species diagnosis, evolution and phylogeny research of other tephritid fruit flies in the future.

The karyotype, genome survey, and assembly of Mud artemisia (Artemisia selengensis).

BACKGROUND: Artemisia selengensis is traditional Chinese medicine and phytochemical analysis indicated that A. selengensis contains essential oils, fatty acids and phenolic acids. The lack of reference genomic information may lead to tardiness in molecular biology research of A. selengensis. METHOD AND RESULTS: Karyotype analysis, genome survey, and genome assembly was employed to acquire information on the genome structure of A. selengensis. The chromosome number is 2n = 2x = 36, karyotype formula is 28 m + 8Sm, karyotype asymmetry coefficient is 58.8%, and karyotypes were symmetric to Stebbins' type 2A. Besides, the flow cytometry findings reported that the mean peak value of fluorescent intensity is 1,170,677, 2C DNA content is 12 pg and the genome size was estimated to be approximately 5.87 Gb. Furthermore, the genome survey generates 341,478,078 clean reads, unfortunately, after K-mer analysis, no significant peak can be observed, the heterozygosity, repetitive rate and genome size was unable to estimated. It is speculated that this phenomenon might be due to the complexity of genome structure. 37,266 contigs are preliminary assembled with Oxford Nanopore Technology (ONT) sequencing, totaling 804 Mb and GC content was 34.08%. The total length is 804,475,881 bp, N50 is 29,624 bp, and the largest contig length is 239,792 bp. CONCLUSION: This study reveals the preliminary information of genome size of A. selengensis. These findings may provide supportive information for sequencing and assembly of whole-genome sequencing and encourage the progress of functional gene discovery, genetic improvement, evolutionary study, and structural studies of A. selengensis.

Meiotic recombination mirrors patterns of germline replication in mice and humans.

Genetic recombination generates novel trait combinations, and understanding how recombination is distributed across the genome is key to modern genetics. The PRDM9 protein defines recombination hotspots; however, megabase-scale recombination patterning is independent of PRDM9. The single round of DNA replication, which precedes recombination in meiosis, may establish these patterns; therefore, we devised an approach to study meiotic replication that includes robust and sensitive mapping of replication origins. We find that meiotic DNA replication is distinct; reduced origin firing slows replication in meiosis, and a distinctive replication pattern in human males underlies the subtelomeric increase in recombination. We detected a robust correlation between replication and both contemporary and historical recombination and found that replication origin density coupled with chromosome size determines the recombination potential of individual chromosomes. Our findings and methods have implications for understanding the mechanisms underlying DNA replication, genetic recombination, and the landscape of mammalian germline variation.

Streptomyces coffeae sp. nov., an endophytic actinomycete isolated from the root of Coffea arabica (L.).

An endophytic actinobacterium, designated strain CA1R205T, was isolated from the surface-sterilized root of Coffea arabica L. collected from Ratchaburi province, Thailand. The taxonomic position of this strain was evaluated using a polyphasic approach. The strain produced light yellowish brown to dark brownish black substrate mycelium and greyish white aerial mycelium. The spiral spore chains were produced directly on aerial mycelium. CA1R205T was found to have ll-diaminopimelic acid in the cell peptidoglycan, galactose, glucose, mannose and ribose as whole-cell reducing sugars, MK-10(H4), MK-9(H6), MK-10(H2), MK-9(H4), MK-10(H6) and MK-10(H8) as menaquinones and iso-C15 : 0, anteiso-C15 : 0, iso-C16 : 0 and C16 : 0 as major fatty acids. Diphosphatidylglycerol, phosphatidylethanolamine, hydroxyphosphatidylethanolamine, phosphatidylglycerol and phosphatidylinositol were detected in the cells. These characteristics were consistent the typical chemotaxonomic properties of members the genus Streptomyces. The taxonomic affiliation at the genus level of this strain could be confirmed using its 16S rRNA gene sequence data. CA1R205T showed the highest 16S rRNA gene sequence similarity value to Streptomyces rapamycinicus NRRL B-5491T (98.9 %), followed by Streptomyces iranensis HM 35T (98.8 %). Digital DNA-DNA hybridization and average nucleotide identity-by blast (ANIb) values between CA1R205T and S. rapamycinicus NRRL B-5491T were 27.2 and 81.5 %, respectively. The DNA G+C content of genomic DNA was 70.7 mol%. Due to the differences in physiological, biochemical and genotypic data, CA1R205T could be discriminated from its closest neighbour. Thus, CA1R205T should be recognized as representing a novel species of the genus Streptomyces, for which the name Streptomyces coffeae sp. nov. is proposed. The type strain is CA1R205T (=TBRC 11244T=NBRC 114295T).

Budding yeast Rap1, but not telomeric DNA, is inhibitory for multiple stages of DNA replication in vitro.

Telomeres are copied and reassembled each cell division cycle through a multistep process called telomere replication. Most telomeric DNA is duplicated semiconservatively during this process, but replication forks frequently pause or stall at telomeres in yeast, mouse and human cells, potentially causing chronic telomere shortening or loss in a single cell cycle. We have investigated the cause of this effect by examining the replication of telomeric templates in vitro. Using a reconstituted assay for eukaryotic DNA replication in which a complete eukaryotic replisome is assembled and activated with purified proteins, we show that budding yeast telomeric DNA is efficiently duplicated in vitro unless the telomere binding protein Rap1 is present. Rap1 acts as a roadblock that prevents replisome progression and leading strand synthesis, but also potently inhibits lagging strand telomere replication behind the fork. Both defects can be mitigated by the Pif1 helicase. Our results suggest that GC-rich sequences do not inhibit DNA replication per se, and that in the absence of accessory factors, telomere binding proteins can inhibit multiple, distinct steps in the replication process.

A human respiratory tract-associated bacterium with an extremely small genome.

Recent advances in culture-independent microbiological analyses have greatly expanded our understanding of the diversity of unculturable microbes. However, human pathogenic bacteria differing significantly from known taxa have rarely been discovered. Here, we present the complete genome sequence of an uncultured bacterium detected in human respiratory tract named IOLA, which was determined by developing a protocol to selectively amplify extremely AT-rich genomes. The IOLA genome is 303,838 bp in size with a 20.7% GC content, making it the smallest and most AT-rich genome among known human-associated bacterial genomes to our best knowledge and comparable to those of insect endosymbionts. While IOLA belongs to order Rickettsiales (mostly intracellular parasites), the gene content suggests an epicellular parasitic lifestyle. Surveillance of clinical samples provides evidence that IOLA can be predominantly detected in patients with respiratory bacterial infections and can persist for at least 15 months in the respiratory tract, suggesting that IOLA is a human respiratory tract-associated bacterium.

Genomic GC content drifts downward in most bacterial genomes.

In every kingdom of life, GC->AT transitions occur more frequently than any other type of mutation due to the spontaneous deamination of cytidine. In eukaryotic genomes, this slow loss of GC base pairs is counteracted by biased gene conversion which increases genomic GC content as part of the recombination process. However, this type of biased gene conversion has not been observed in bacterial genomes, so we hypothesized that GC->AT transitions cause a reduction of genomic GC content in prokaryotic genomes on an evolutionary time scale. To test this hypothesis, we used a phylogenetic approach to analyze triplets of closely related genomes representing a wide range of the bacterial kingdom. The resulting data indicate that genomic GC content is drifting downward in bacterial genomes where GC base pairs comprise 40% or more of the total genome. In contrast, genomes containing less than 40% GC base pairs have fewer opportunities for GC->AT transitions to occur so genomic GC content is relatively stable or actually increasing. It should be noted that this observed change in genomic GC content is the net change in shared parts of the genome and does not apply to parts of the genome that have been lost or acquired since the genomes being compared shared common ancestor. However, a more detailed analysis of two Caulobacter genomes revealed that the acquisition of mobile elements by the two genomes actually reduced the total genomic GC content as well.